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Abstract 

Anomalous dimensions of high- twist Wilson operators have a nontrivial scaling behavior in the 
limit when their Lorentz spin grows exponentially with the twist. To describe the corresponding 
scaling function in planar Af = 4 SYM theory, we analyze an integral equation recently proposed 
by Freyhult, Rej and Staudacher and argue that at strong coupling it can be casted into a form 
identical to the thermodynamical Bethe Ansatz equations for the nonlinear 0(6) sigma model. 
This result is in a perfect agreement with the proposal put forward by Alday and Maldacena 
within the dual string description, that the scaling function should coincide with the energy 
density of the nonlinear 0(6) sigma model embedded into AdSs x S^. 
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1 Introduction 



Anomalous dimensions of high-twist Wilson operators have a nontrivial scaling behavior in the 
limit when their Lorentz spin grows exponentially with the twist Recently, Alday and Mal- 
dacena [2] used a dual stringy description of such operators in A/" = 4 super- Yang-Mills theory 
(SYM) to put forward the proposal that, at strong coupling, the corresponding scaling function 
should coincide in a suitable limit with the energy density of a two-dimensional bosonic 0(6) 
sigma-model. In this paper, we establish the same relation on the gauge theory side of the 
AdS / CFT correspondence. 

The operators under consideration are built from L complex scalar fields X(0) and N covariant 
derivatives = {n ■ D) projected onto the light-cone direction = 



with N = ki + . . . + kL- These operators mix under renormalization and their anomalous 
dimensions, defined as eigenvalues of the mixing matrix, have a rich structure. For given N and 
L, they occupy a band and the properties of anomalous dimensions are different in the lower 
and upper part of the spectrum [3]. In what follows we shall study the minimal anomalous 
dimension 'jN,L{g) in the planar A/" = 4 SYM theory, both at weak and strong coupling, in the 
limit N,L oo. The motivation for considering this particular limit is twofold and it goes 
beyond the scope of A/" = 4 model. High spin/twist Wilson operators analogous to ([T]) naturally 
appear in QCD (with a different partonic content though) in the operator product expansion 
description of deeply inelastic scattering in the semi-inclusive regime, the so-called large Bjorken 
X limit [1]. On the other hand, in planar J\f = 4 SYM, the Wilson operators ([1]) admit a dual 
description in terms of folded strings spinning on the AdSs x S^ [5], [6] and their anomalous 
dimensions can be found at strong coupling from the AdS/CFT correspondence [7]. 

The minimal anomalous dimension 'yN,L{g) is a complicated function of the Lorentz spin A^, 
twist L and 't Hooft coupling g"^ = A/(47r)^ = (7YM^c/(47r)^. For finite twist L and large spin A^ 
it has a remarkable logarithmic (Sudakov) scaling behavior [S", "3] 



where the coupling dependent prefactor does not depend on the twist L and is given by the 
cusp anomalous dimension [9l |8]. The situation changes however when A^ and L become large 
simultaneously. In that case, the anomalous dimension '^N,L{g) still scales logarithmically for 
sufficiently large A^. However, in distinction with ([2]), the value of the spin for which transition 
into the Sudakov regime takes place now depends both on the twist L and on the coupling 
constant. As was shown in Ref. [1], different regimes of the anomalous dimension 'jN,L{g) at 
large L and A^ ^ L are controlled by the order parameter ^, which is given at weak coupling by 
^{g <1) = \yi{N/L)/L and at strong coupling by ^(^f > 1) = g\n{N/L)/L: 

• For ^ < 1 the minimal anomalous dimension has a BMN like [10] scaling both at weak 
and strong coupling [6l [121 [I] 



(1) 



jNAg) = ^^cnsp{g) lnN + 0{N') 



N oo, L = fixed , 



(2) 



7iv,L = 8g' \n\N/L)/L + 0{gHn\N / L) / L^) 



A^, L -> oo , C < 1 • 



(3) 



^As was shown in Ref. [Tl], this scahng breaks down at higher orders in the BMN coupling jl? . 
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• For ^ 1 the minimal anomalous dimension grows logarithmically with and it has the 
following scaling behavior both at weak and at strong coupling [H [TT|, [131 El 

-fN Aa) = I'^^cnspig) + e{g, j)]\n N + ... , AT, L ^ oo , j = ^ = fixed . (4) 

Here e{g,j) is a nontrivial function of the scaling variable j satisfying e{g,j = 0) = and ellipses 
denote terms suppressed by powers of 1/L. We note that the condition C, ^ 1 automatically 
implies that j <^ 1 at weak coupling and j <^ g aX strong coupling. 

At weak coupling, the scaling function e{g,j) can be found in a generic (supersymmetric) 
Yang-Mills theory in the planar limit by making use of the remarkable property of integrability 
[151 [16]. Ill QCD, the dilatation operator in the SL{2) sector is integrable for the special class of 
maximal helicity operators to two loops at least [T71 [TH [19] , whereas in A/" = 4 SYM theory it is 
believed that integrability gets extended in this sector to all loops [20], [21]. The scaling function 
has the following form at weak coupling 

<g.3) = eiig)j + ^2ig)f + e3{g)f + • • • , (5) 

with the coefficient functions ek{g) given by series in g^. To one-loop order, one finds in planar 
J\f = A theory [I] | 

ei = 4 [As) - A^s)] g' + ..., = ■ g' + . . . , e, = ^(-Ais))g' + ... (6) 

where s = 1/2 is the conformal spin of scalar fields entering ([1]) and ip^x) = (lnr(a;))' is the Euler 
psi-function. Recently, an integral equation has been proposed by Freyhult, Rej and Staudacher 
(FRS) [13] to describe the function e{g,j) for arbitrary values of the scaling parameter j and the 
coupling constant g. It generalizes the BES equation [22] which governs the dependence of the 
cusp anomalous dimension Tcuspig) on the coupling constant in planar Af = 4 SYM. One of the 
consequences of this equation is that 62(5') = to any order in g"^. 

At strong coupling, the gauge/string correspondence relates the minimal anomalous dimension 
to the energy of a folded string spinning on the AdSs x background [Sj, [B] 

A = N + L + ^mA9), (7) 

with A^ and L being angular momenta on AdSa and S^, respectively. Semiclassical quantization 
of this state yields the expansion of the scaling function e{g,j) in powers of 1/g. The first 
two terms of this expansion have been computed in Refs. [231 [ISl [21]. Recently, Alday and 
Maldacena [2] put forward the proposal that the scaling function e{g,j) can be found exactly at 
strong coupling in the limit j <^ g and j/m = fixed (with the parameter m defined below in Q). 
They argued that quantum corrections in the AdSs x sigma model are dominated in this limit 
by the contribution of massless excitations on whose dynamics is described by a (noncritical) 
two-dimensional bosonic 0(6) sigma-model equipped with a UV cut-off determined by the mass 
of massive excitations. 



^Our definition of e{g,j) and j is slightly different compared to that of Ref. [2]. 

•^We displayed in ^ the dependence on s to indicate that the same expressions hold in a generic (planar) 
Yang-Mills theory in integrable SL{2) sector of aligned helicity operators built from gauge fields (s = 3/2) and 
gaugino (s = 1). 
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The 0(6) sigma model emerges within the AdS/CFT as an effective two-dimensional theory 
describing the scaling function e{g,j) in the AdS5 x sigma model. More precisely, the quantity 



eo(6) = ^ = Jim - T^^M (8) 

has the meaning of the energy density in the ground state of the 0(6) model corresponding 
to the particle density p = L/{2\nN) = j/2. The 0(6) sigma model is an exactly solvable 
theory [23]. It has a nontrivial dynamics in the infrared and massless excitations acquire mass 
through dimensional transmutation mechanism [26l [27]. The exact value of the mass gap was 
found using the Bethe Ansatz in Ref. [28]. In terms of parameters of the underlying AdSs x 
sigma model, its value reads [2] 

m = kg'/' e-^ [1 + 0(1/ g)] , k = . (9) 

The dependence of m on the coupling constant is fixed by the two-loop beta-function of the 
0(6) model whereas the prefactor k was determined in [2] by matching the first few terms of 1/g 
expansion of eo{6) into semiclassical expansion of €{g,j) computed in [TT] . 

The appearance of mass gap in the dual stringy description of the scaling function has dra- 
matic consequences for A/" = 4 SYM at strong coupling. It suggests that the strong coupling 
expansion of the anomalous dimensions in A/" = 4 SYM should receive nonperturbative correc- 
tions characterized by a new 'hidden' scale m. This is a rather unexpected and surprising feature 
given the fact that the anomalous dimensions capture dynamics at short distances and their weak 
coupling expansion is free of nonperturbative corrections. Indeed, such nonperturbative correc- 
tions have been identified in Ref. ^29j in the strong coupling expansion of the cusp anomalous 
dimension 



oo 



rcusp(^) = 2^ Ck/g" + am' + o{m') , (10) 

k=-l 

where the perturbative coefficients Ck grow factorially at large k and the value of a depends 
on regularization of Borel singularities in the perturbative series. The 0{m^) corrections are 
exponentially small at strong coupling but they provide a leading contribution in the region 
5f ~ 1 in which the transition from the strong to weak coupling regime takes place. 

In this paper, we solve the FRS equation [l^ using the approach developed in Ref. [29] and 
evaluate the scaling function e{g,j) in A/" = 4 SYM theory at strong coupling. More precisely, we 
demonstrate that, in the limit g ^ oo and j/m = fixed, the integral equation for e{g,j) can be 
casted into a form identical to the thermodynamical Bethe Ansatz equations for the nonlinear 
0(6) sigma model [2B]. This allows us to determine e{g,j) at strong coupling for different values 
of the ratio j/m and to translate nontrivial properties of the 0(6) sigma model (mass gap 
generation at large distances and asymptotic freedom at short distances) into the corresponding 
scaling behavior of e{g,j): 

• For j <^ m <^ g, or equivalently A^ ^ e^/'*", 

ij,9)+J=m' - + + 0{jVm') . (11) 

m 24 \m J ^ ' 
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In this expansion, 0{j) terms are in agreement with the numerical solution to the FRS 
equation found in [30]. Also, O(j^) term is absent, in a striking similarity with vanishing 
of the function €2(5') in the weak couphng expansion and (jSD- Note that the nonpertur- 
bative corrections in ffTOj) and ffTTl) are of the same order in and are added together in 
the expression for the anomalous dimension (jl]). 



For m <^ j <ti g, or equivalently e^^"^ ^ N ^ e^^^, 



9)+j= J 



, ^/'lnln(j7m) 



;ln(j/m) \ In {j/m 



(12) 



This expression resums through renormalization group (an infinite number of) perturbative 
corrections inl/g with coefficients proportional to and enhanced by powers of ln{j/g) [H 



• For j <^ g and j/m = fixed, or equivalently N ~ e^/*", the function e(j, g) does not admit 
a simple representation and it can be found as the solution to thermodynamical Bethe 
Ansatz equations for the 0(6) model |28] (see Eqs. fl65l) - flTOj) below). 

These results are in a perfect agreement with the Alday-Maldacena proposal [2] and, therefore, 
constitute a nontrivial test of the AdS/CFT correspondence. 

The paper is organized as follows. In Sect. 2, we reformulate the integral equation proposed 
in Ref. fL^ and argue that it can be significantly simplified by employing a nontrivial change of 
variables found in [25]. In Sect. 3 we work out the expansion of the scaling function at small j. We 
demonstrate that at weak coupling it agrees with the known one- loop result, whereas at strong 
coupling it matches the string theory prediction by Alday and Maldacena including expression 
for the mass gap Q. In Sect. 4 we show that the integral equation for the scaling function can 
be mapped for g ^ 00 and j/m = fixed into thermodynamical Bethe Ansatz equations for the 
energy density in the ground state of two-dimensional 0(6) sigma-model with the particle density 
j/2. Section 5 contains concluding remarks. Some details of our calculations are summarized in 
two appendices. 



2 The scaling function in AT = 4 SYM 

The derivation of the scaling function relies on integrability of the dilatation operator in the 
SL{2) sector. To one-loop order, the function e{g,j) was determined in Ref. pLj from the analysis 
of Bethe Ansatz equation in the scaling limit N, L —>■ oo and j = fixed. The Bethe Ansatz 
solution for e{g,j) is characterized in this limit by two sets of parameters, the Bethe roots and 
the so-called small roots of the transfer matrix, which describe excitations dubbed 'magnons' 
and 'holes', respectively. As was shown in [1], both sets of parameters form a dense distribution 
on the real axis with the holes confined to the interval [—a, a] and magnons to the union of two 
intervals [—00,— a] U [a, cx)]. The scaling function is uniquely determined by the corresponding 
distribution densities of holes and magnons. 

This analysis was recently extended to all loops in [13] by employing the asymptotic Bethe 
Ansatz approach proposed in Refs. [311 |2T]. It led to the FRS equation 

a{t) = -^(^K:{t,0) -4 dt' t{t,t')a{t')^ , (13) 
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whose solution a^t) = a{t] g,j) is related to the scaling function as 

f{g,j) = 2r,,3p(^?) + e{g,j) = j + 16a(0) . (14) 

The kernel ]C{t,t') is given by a rather complicated expression that can be found in Ref. [H] 
(see also Eq. ( l22l) below). It takes into account a nontrivial scattering phase which satisfies the 
crossing symmetry [32] and whose explicit form was found in [33] . Later in the paper we shall use 
another, equivalent formulation of the integral equation (1T3!) which is more suitable for studying 
the strong coupling limit. 

For j ^ the integral equation (|T3l) reduces to the BES equation [22] and its solution 
a{t]g,j = 0) determines the cusp anomalous dimension f{g,j = 0) = 2rcusp(5')- Then, the 
scaling function e{g,j) admits the representation 

<9,j) = f{9,j)-f{9,0)- (15) 

At weak coupling the function f{g,j) admits an expansion in powers of g"^ and it vanishes as 
(7^ — > 0. Then, it follows from f|T^ that the integral equation f|T3l) has a nontrivial solution at 
g = satisfying a{0; g = 0, j) = -j/16. 



2.1 Integral equation 

It is convenient to split a{t) into a sum of three terms 



a(t) 



e-*/'7h(t) + ^7(2^?t)-^Jo(2^?t) 



(16) 



with Jo{x) being the Bessel function, and rewrite (fT3|) as a system of coupled integral equations 
for the functions 7h(^) and 'j{2gt) 



POO 

7h(t) = Kh(t, 0) - 4 / dt' Mt, t') e*'/2 a{t') 
Jo 

2t I f°° 
^{2gt) = — K(t,0)-4 dt' K(t,t')a(t') , 
9 1 Jo 



(17) 



with the kernels K-^ and K defined below in Eqs. (ITSl) and (fT9l) . respectively. The rationale 
behind the decomposition (fT6l) is that, as we will show below, the functions 7h(^) and 'y{2gt) 
have different analytical properties. Namely, the Fourier transforms of 7h(^) and 'y{2gt) defined 
as solutions to (ITTj) have the support on the intervals [—a, a] and [—2g, 2g], respectively, with the 
parameter a depending on g and j. 

The first relation in f|T7|) involves the kernel 

t cos {at') sin (at) — t' cos (at) sin {at') , , 



Mt,t') 



27r(t2 - 
sin(a(t-t')) ^ sm{a{t + t')) 



1 

47r 



t-t' t + t' 

and the kernel in the second relation in ffTTI) is defined as 



K{t,t')=g'K^'\2gt,2gt')-Ag' / dt" K^'\2gt,2gt")K^{t" ^t')^^''-'"^^ (19) 
K^'\t,t') = K4t,t')+K.{t,t')+8g' / 4^ir_(t,20ir+(2<^t",t'), 
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where K^'^\t,t') coincides with the BES kernel [22] and parity even/odd kernels K±(~t,t') = 
K±{t, -t') = ±K±{t,t') are given by 

^ tJ,(t)Jo(t2-t;Jo(t)Jx(tO ^ 2 ^^2, _ l)J,„_,(t)J,„_,(0 , (20) 



t2 - t'2 

t'Ji{t)Mt')-tMt)J,{t') 2 



n>l 



n>l 



Finally, the parameter a = a{g,j) is related to the solution of the integral equation as 

4a 16 f°° ^^^^ sin(at) 



J = dta{t)e'""-^. (21) 

Being combined together, the relations f|T6|) - fl2T]) are equivalent to the integral equation f|T3l) 
with the kernel given by 

^(t, t') = K{t, t') + t-iKh(t, t') e(*'-*)/2 _ ^o(2^t) sin K) ^.72 . (22) 

We would like to stress that the functions 7(t) and "^^lit) depend on the coupling constant g and 
the scaling variable j. The relations fll7p - 0211) simplify significantly as j — 0. In this limit, 
we have a ~ j so that the kernel -ft'h(^) ^') and the function 7h(^) vanish simultaneously as j ^ 
leading to }Cj=o{t,t') = g'^K^^\2gt,2gt') and to the BES equation for the function '-fj=o{2gt). 

To determine the scaling function f|T^ . we have to solve the integral equations f|T7|) and, then, 
substitute the solution for cr{t), Eq. f|T6l) . into f|T^ . Let us apply f|T6l) to evaluate o-(O). We put 
t = in both sides of the first relation in f[T7|) and take into account the relations flTSl) and fl?Il) 
to verify that 7h(0) = j/S for arbitrary j and g. Substituting this relation into ffTB]) and ffn|) we 
find the following representation for the scaling function 

2.2 Parity decomposition 

Similar to the BES equation [Ml [35]) the integral equations (JTTj) can be further simplified by 
separating expressions on both sides of f|T7|) into terms with a definite parity under t — > — t. 

According to definition f|T8|) . the kernel K]^{t,t') is an even function of t and t'. Then, it 
immediately follows from the first relation in f|T71) that 7h(t) is an even function of t 

7h(-t) = 7h(t) , 7h(0)=j78, (24) 

where the second relation holds for arbitrary coupling g. As follows from (fT9l) . the kernel K{t, t') 
does not have a definite parity under t ^ ~t and, as a consequence, the function 7(t) can be 
decomposed as 

7(t) = 7+(t) + , 7±(-t) = ±7±(t) ■ (25) 
We substitute fl25l) into the second relation in f|T7|) and match parity even/odd terms in both 
sides of f|T7j) to obtain after some algebra the system of coupled equations for the functions 7±'''^^ 

7T(2^t) 



2gt 



POO 

-4/ dt'K±{2gt,2gt') e-''/^^^{t') + ait') +b±{2gt), (26) 

Jo L -I 



"^In comparison with the second relation in (|17p . these equations do not involve the convolution of the kernels 
K'-°^ and Kh which appears in the kernel K(t, t'), Eq. (|19p . 
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where the kernels K± are defined in (l20l) . a{t) is given by ( fT6l) and the notation was introduced 
for 

b4t)=K^it,0), b4t)=Ag 



e*' -1 



K4t,2gt')-f42gt') . 



(27) 



We recall that the kernels K± are given by sum over the product of Bessel functions, Eq. (!20|) . 
Substitution of (120!) into (!26l) yields the expression for 7±(t) in the form of Neumann series over 
Bessel functions 



7_(t) =2^ (2n-l)J2n-i(t)7: 

7+(t) =2^ (2n) J2n(t)72n 



2n-l , 



(2^ 



with the expansion coefficients 7^, = 'jn{g,j) given by the following expressions 



7n 



If 



dt Jn{2gt) 



t 



[7+(2^?t)-(-iyV(2^7t)]. 



Here the notation was introduced for the coefficients /i„ 

_2 r^dt U2gt) 
9 Jo t e* -1 



e'/'j^{t)-iM2gt) 



(29) 



(30) 



We recall that the scaling function fl23|) is determined by the asymptotic behavior of the function 
7(t) for t ^ 0. Substituting ( !25l) and f l28l) into ( l23l) and taking into account that Jn(t) ~ at 
small t, we find the scaling functions f |T4l) and ( |T5l) 



/(^,j) = 1%V(^,J; 



e((7,j) = i6^?' biig,j)-iiig,o)] 



(31) 



where 71 is the lowest order coefficient in the series ( |28l) . 7-(t) = li{g,3)t + O(t^). 

We can obtain another representation for the expansion coefficients 7„ in fl28l) by making use 
of orthogonality of the Bessel functions 



/"°° dt [°° dt 1 

(2n - 1) / — J2n-l{t)J2m-l{t) = 2n — J2n{t)J2m{t) = -S„ 

Jot Jo 't ^ 



In this way, we obtain from fl28l) 

dt 



l2n-\ 



t 



J2n^l{t)-f-{t) , 72„ 



dt 



J2n{t)l+{t) 



(32) 



(33) 



Here we have tacitly assumed that the Bessel series on the right-hand side of (!28|) are convergent 
on the real axis and that the sum over n can be interchanged with the integral over Matching 
(133|) into (!29|) we obtain an (infinite) system of equations for the functions 7±(t) 







°° dt 

T 
~ dt 
T 



J2n-l{'2gt) 



7_(2^t) 7+(2^t) 
1 - e-* e* -1 



J2n(2f7t) 



7+(2(?t) 7-(2^7t) 



2 

h2nig,j 



Sn,l + h2n~l{g,j) 



(34) 



^In particular, performing the Fourier transformation on both sides of (j28p and using weU-known properties 
of the Bessel functions, we find that the Fourier transforms of the functions ^±{2gt) have support on the interval 
[-25,25]. 
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valid for n ^ 1. These relations should be supplemented with expressions (1301) for inhomogeneous 
terms hn- The latter depend on the function 7h(^) which depends on its turn on 7(t) = 7+(t) + 
7_(t) and satisfies the integral equation (IT7|1 



7h(t) =i^h(t,0)-2 



dt'Khit,t') 
sinh(t'/2) 



U(2gt')+e-''/'^^it')-lM2gt') 



(35) 



with the kernel Kh(t,t') given by f|T8|) . 

Let us now examine the system of integral equations flM|) and fl5^ for j = 0. In this case, the 
kernel Kh{t,t') vanishes and we deduce from fl5S]) and fl5U]) that 7h(i; j = 0) = hn{g,j = 0) = 0. 
Then, the relations (!34|) are reduced to the BES equation [22] and we denote their solution as 



7±(t; J = 0) ^ 7^ (t) =tr,nsA9)/M + 0(t' 



(36) 



The functions 7^°'* (t) were constructed at weak coupling in Refs. [221 EE] and at strong coupling 
in Refs. [33 [371 [291 [381 [39] . 

To evaluate the scaling function f[3T]) for arbitrary j, it is sufficient to determine the coefficient 
71 defined in (1281) rather than the function 7-(t) itself. This can be achieved by using the fact 
that the solutions to the integral equations flMl) and f[35l) satisfy Wronskian-like relations. These 
relations allow us to obtain the following relation for the scaling function (see Appendix A for 
details) 



di 7f (2^t)-7i°H2^t) 



J/2 



7h(t; j) 



4,, r^7f(2,t). 
Jo >^ 



(37) 



It is remarkable that e{g,j) only depends on the solution of the BES equation, 'y^\2gt), and 
the 'hole' function 7h(^,j)- It is the latter function that generates a nontrivial dependence of 
^id^j) oil the scaling parameter j. We note that the integrals on the right-hand side of ( 1371) are 
well-defined for t in virtue of 7!°^) ~ t, -ff\t) ~ and 7h(0; j) = j/8, Eqs. ([281) and ([Ml), 
respectively. 



2.3 Change of variables 



The system of integral equations (1341) has been analyzed in Refs. [291 [39] in the special case /i„ = 
corresponding to j = 0. The starting point in [29J was the change of variables 7±(t) — > r±(t) 
that eliminated the dependence of the BES kernel on the coupling constant. It is remarkable 
that the same change of variables also applies to ( l34l) for arbitrary j 



r±(t; j) = 7±{t-j) T coth ( ^ ]i^{t;j) 



or conversely 



27±(t;j) 



r±(t;j)±tanh( ;^VT(i;j) 



cosh(t/2^)y \2g J 
Upon this change of variables, the relations ([Ml) take the form (for n > 1 



dt 



Jn{t) 



r_(t) + (-i)"r+(t) 



5nA + 2hr^ 



(38) 



(39) 



(40) 
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where /;,„— coefficients are given by (!30l) and the integral kernel on the left-hand side is the same 
as for j = 0. It is convenient to apply the Jacobi- Anger expansioii^ (with u = sinip and (f real) 



n>l 



{ J2„(t) [cos(2nv9) - 1] + iJ2n-i(t) sin ((2n - l)(p)} 



(41) 



and replace an infinite system of equations ( l40l) into a single relation depending on the m— parameter 

r-oo ^2igtu 



dt 







-itu 



2 - 



e*-l 



(42) 



It is important to stress that this relation only holds for — 1 < m < 1. Here r-|-(t) and 7h(i^) are 
real functions of t. Taking real and imaginary part on both sides of fH2l) . we can replace P2l) by 
a system of two coupled equations. Later in the paper we shall make use of this system. 

For j = the integral on the right-hand side of ( 142|) vanishes. Its solution at strong coupling, 
r±(t; j = 0) = r|t°^(t), was constructed in Refs. [291 [39]. For j ^ the relation ([42D should be 
supplemented with integral equation for the function 7h(t). 



3 Scaling function at small j 

To find the scaling function ( 1371) . we have to solve a complicated system of coupled integral 
equations for the functions 7h(^) and 7±(t). In this section we shall evaluate the scaling function 
e{g,j) at small j and generalize consideration to arbitrary j in the next section. As we will see 
in a moment, the ffist few terms of the small j expansion of e{g,j) can be found exactly without 
solving the integral equations for the functions 7±(t). 



3.1 Small j expansion 



To begin with let us consider the integral equation (|T7j) for the function 7h(i)- It involves the 
kernel K\i{t,t') defined in (JTSl) . According to its definition, KY,{t,t') depends on the parameter a 
which depends on its turn on j and satisfies the relation fl2Tl) 



4a 8 dt sin(at) 



TT 



IT 



t sinh(t/2) 



3 



(43) 



We recall that the parameter a determines the interval [—a, a] on which holes rapidities condense 
in the limit N,L oo and j = fixed. As was shown in [T], for j — > this interval shrinks into a 
point a ~ j and the scaling function vanishes e{g,j = 0) = 0. 

Assuming that a ~ j we expand the kernel K^^lt, t') in powers of a and obtain from (|T8l) 

sin(at') 



1 
2^ 



6 



2711' 



0(a' 



(44) 



where expansion runs in odd powers of a only. We substitute this relation into the right-hand side 
of (135|) and observe that the t'— integral there coincides with analogous integral on the right-hand 
side of (1431) . Combining the two relations together, we find from ( |35l) 



7h(t) = i 



l--{atf + 0{a') 



(45) 



'More precisely, we differentiate botli sides of this relation with respect to u and, then, divide by t. 
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in agreement with (!24l) . 

To determine a we return to the relation ( H3il and expand its both sides in powers of a and 
j. We notice that the last two terms in the square brackets there vanish linearly as j — ^ while 
the first term is given by 'j{2gt) = 'y^^\2gt) + 0{j). In this way, we find from fH3|) 



a = ^ + 0{f) , (46) 



where the (7— dependence resides in the normalization factor 



. = 2-2,/ *2f(M±^, (47) 

^ ' sinh(t/2) ^ ^ 







Substitution of (H6|l into (H5|) yields the expansion of 7h(i^) at small j. 

According to (!37j) . the j— dependence of the scaling function e{g,j) is controlled by the func- 
tion 7h(t). Replacing 7h(^) in (j37l) by its small j expansion (l45l) . we get 

<9,j)=jei{g)+fesig) + 0{f), (48) 

where the coefficient in front of equals zero for any g and the coefficient functions £1(5') and 
es^g) are given by 



e.(.)^-4.y^ - 



lf\2gt) j^^\2gt) 



eM--4^ r,,,7f(2,t)-7i^2,t) 



(49) 



12/«2 7o sinh(t/2) 
We would like to stress that these relations hold for arbitrary coupling g. 

3.2 Weak coupling 

We recall that the functions 'y±\t) satisfy the BES equation. At weak coupling, this equation 
can be solved by iterations leading to [22|, [36] 



7^-\t) = (l - Mt) + O {g') , (50) 

7f (t)=4/C3J2(t)+0(/). 
Plugging these expressions into (H9l) and expanding €i{g) and e^i^g) in powers of g"^ we get 

ei (g) = -Sg' In 2 + / Qtt^ In 2 + I6C3) + 0(^7^) , (51) 

es{9) = /^^^Cs + 9' (I As - fvr^Cs) + 0(^7^) . 

It is easy to see that the one-loop corrections to ei{g) and £3(5') are in an agreement with the 
relation for s = 1/2. 
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It is straightforward to evaluate higher order corrections to ei{g) and e^lg) by taking into 
account subleading terms on the right-hand side of (1501) . It is interesting to notice that higher 
order corrections to ti{g) proportional to 'In 2' are controlled to all loops by the cusp anomalous 
dimension 

ei((7) = -2reusp(^7)ln2 + ..., (52) 

Indeed, such term originates from the second integral on the right-hand side of ( l49l) with 7!°'' {2gt) 
replaced by its leading small— t asymptotic behavior ( l36i) . 

The relations (15T!) and (152!) are in agreement with the results of [H]. 



3.3 Strong coupling 

To evaluate ei{g) and £3(5') at strong coupling, it is convenient to substitute 'y±\t) on the right- 
hand side of fH9l) by their expressions fl39|) in terms of the functions r|^'*(t) = T±{t;j = 0). 
This leads to the following representation for ei{g) 



ei{g) = -2g 



dt 
T 



cosh 

cosh(t/2^) J 



cosh {t/2g) \ ' ^ ' 



rf(t) 



(53) 



The functions T±{t;j) satisfy the integral equation ( H2i) . For j = the integral on the right-hand 
side of ( l42i) vanishes and the resulting equations for rl^-* (t) can be written as 



dt sm{ut) rL°^(t) + rf (t) 



0, 



(54) 



dt cos{ut) rL°^(t) -rV''(t) 



where m— parameter satisfies the condition u"^ < 1. 

We observe that the relations (!53|) involve the same combinations of r|^^ (t) as in (15^ . To make 
use of 1^^, it is suggestive to replace the two factors on the right-hand side of (153|) involving the 
ratios of hyperbolic functions by their Fourier integrals. Going through calculation (the details 
can be found in Appendix B) we find 



e,{g) = -AV2g' 



sinh [gnu) 



+ 



+ 



du u 

_i cosh [zgnu) 

°° cosh (gnu) 
du 



(55) 



du 



cosh {2g7ru) 
sinh {gnu) 
cosh {2gnu) 







°° dt . , , 
— sm [ut) 



cos [ut) 







t 



rL°^(t) -r 



rL°)(t) + rf (t) 



Here we split the Fourier integral into sum of two terms corresponding to < m < 1 and 
1 < M < 00 and applied ( IMl) in the first term. 

Let us now examine the relation (155!) at strong coupling. It is straightforward to work out the 
asymptotic expansion of the first term on the right-hand side of (155|) at large g. In the remaining 
two terms, we replace the ratios of hyperbolic functions by their large g expansion and find after 
some algebra (see details in Appendix B) 



e^{g) = -1 + m + O (e-3f-) 



(56) 
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where the parameter m = m{g) is defined as 



m 



8v^ __8^? 



-TTQ 



Re 



TT 



oo gi(i-^/4) 



(rf (t) + ^rL°^(t) 



We recall that the functions r|^''(t) are solutions to the integral equations (E 



(57) 



According to (!56|) . the function ei{g) does not receive perturbative corrections ml/ g and the 
leading nontrivial correction is given by m which is exponentially small in g. This property has 
the same origin as the appearance of the gap in the distribution of magnon rapidities at strong 
coupling described by the BES equation [2^ ESI ES]- Namely, the distribution density of Bethe 
roots vanishes on the interval [—2g,2g] to each order in perturbative 1/g expansion. These are 
the nonperturbative (exponentially small in the coupling) corrections that generate a nonzero 
distribution density on this interval. 

At large g the integral in fl571) receives a dominant contribution from large t ^ g. In order to 



evaluate (l57Il it suffices to take t = ngr and to replace the functions r^"* (Trgr) by their asymptotic 
behavior for g__-->- oo and r = fixed. The leading asymptotic behavior of T^^^ {t^Qt) in this limit 
reads 



(vr^r)+^rL^^(7r^r 



(0), 



v/2r(3/4)r(l + ^r/4) 
" TT r(3/4 + ir/4) 



du 



u 



u 



1/4 



-iirguT 



+ ... 



(5^ 



where ellipses denote terms suppressed by powers of 1/g. From the point of view of the underlying 
BES equation |39] , this expression refiects nontrivial properties of the scattering phase of magnons 
in planar A/" = 4 S YM theory in the the near fiat space limit [ID] . 

We substitute fl58l) into fl57j) . expand the resulting twofold integral at large g and find after 
some algebra (see Appendix B) 



m 



93/4 1/4 



r(5/4) 



(59) 



This relation is in a remarkable agreement with the expression for the mass gap iQ found by 
Alday and Maldacena [2] from string theory considerations. Moreover, the obtained expression 
for €i{g), Eq. fl56|) . is in agreement with the numerical solution to the FRS equation and its 
analytical estimate [SUl SI] • 

Let us now evaluate e^^g). We start with the second relation in (jlSD and change variables 
following (139|) 



dtt 



^mh^/ (0) 
cosh {t/2g) \ - ^ ' 



ASgK^ 

and analogously for k defined in ( HTll 



rf(t) 



cosh(t/4(7) / (0) (0) 
cosh {t/2g) \}-^^) + ^+ 



(60) 



dt 



cosh jt/Ag) f^m^^. 
cosh {t/2g) \ ' ^ ' 



rf(t) 



(61) 



^To obtain this relation from perturbative solution for F^-* (t) + zF^'' {t) — X]fe>o 9 ''^k{t), one has to take into 
account that Tk{t) ~ i'^Fo(i) at large t and, then, resum an entire series in the double scaling limit g,t ^ 00 and 
t/g — fixed. 
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The integrals entering (!60l) and ( !6T|) look similar to that for ei{g), Eq. ( |53l) . and their calculation 
goes along the same lines as before. Namely, we replace the factors involving the ratios of hyper- 
bolic functions on the right-hand side of fl60l) and (I^Tj) by their Fourier integrals (see Eqs. (IB.ip 
and flB.2|) ). apply the relations fl54l) to evaluate the contribution from the region < u < 1 and 
express the remaining integral over 1 < m < oo in terms of r|^^ (t) . Going through calculation we 
find that the leading contribution to and k as (? ^ oo is proportional to the scale m, Eqs. (I5T 
(see also (1B.5I) ). 



96^2 



-m 



+ O {e-^^s) 



IT 



-m + (e-3"f ) 



(62) 



so that the relation P6l) takes the form a = j /m + . . .. Finally, we combine together the relations 
l8|) . (!56l) and (162!) and obtain small j expansion of the scaling function at strong coupling as 

^2 



e(fi',j) = +m) + 



24m 



(63) 



with m given by flS^ . This relation agrees with the string theory prediction by Alday and 
Maldacena It can be also written as 



+ _( A 

m 24 \ m 



+ 



(64) 



wherefrom we expect that the expansion of the scaling function runs in powers of j/m, or equiv- 



alent ly, efc 



m 



2-k 



at strong coupling. 



4 Scaling function and nonlinear 0(6) sigma model 

For j ~ m, the small j expansion employed in the previous section is not applicable. In this 
section, we will show that for j g and j/m = fixed, the scaling function €{g,j) coincides with 
the energy density of the ground state of the two-dimensional 0(6) sigma model. 

4.1 Exact solution of the 0(6) sigma model 

The exact solution for the ground state energy in the two-dimensional 0(6) sigma-model was 
constructed in Refs. [2S]. It can be summarized as follows. The energy density eo(6) in the ground 
state and the particle density p are given by 

eo(6) = TT / ^^X{0) cosh e, P=7r dOxie), (65) 

where the differential rapidity distribution x(6') has the support on the interval [— -B, B] and 
satisfies the integral equation 

X{e)= ! d9' K{9-e')xi9')+ m cosh e. (66) 

J-B 

Here the kernel K{9) = (In 5'(0))'/(27ri) is related to a logarithmic derivative of the exact 
5*— matrix of the 0(6) model ^25j 



K{9) 



47r2 



271 J V 27r/ V2 2n J \2 2n J cosh^ 



(67) 



13 



where ijj{x) = (lnr(a;))' is the Euler psi-function. Later in this section we will encounter its 
Fourier transform 

2 e* +1 

*<=os(2«t/.)^ 

For p/m ^ 1, or equivalently B —>■ 0, the ground state energy of the model is given by 



K{9) = — I dt cos {2et/7i)-^^^ . (68) 



eo(6) =mp + £-p3 + . . . (69) 

and it coincides with the energy density of a dilute, nonrelativistic Fermi gas with the particle 
density p [28]. 

Comparing flMl) with (1691) we observe that the scaling function at strong coupling is related 
at small j <^ m to the ground state energy of the 0(6) model as[^ 

eo{6) = ^ , ^ = 2 ■ ^ ^ 



We will show in this section that the relation fjTOj) holds for arbitrary j in the scaling limit g ^ oo 
and j/m = fixed. More precisely, we will demonstrate that the FRS equation for the scaling 
function €{g,j) can be written in the form of fl65|l and fl66|) upon identification 



X{0) = - I dt cos {2et/7i)^i,{t), (71) 

or conversely 



7h(t) = ^/ decos{2etMx{e), (72) 

with B = aiT/2. 

4.2 Rapidity distribution 

Let us first demonstrate that the Fourier transform of the function 7h(i) fulfills the same integral 
equation (166!) as the rapidity distribution density for the 0(6) model. 

We recall that 7h(^) satisfies the integral equation (!35|) with the kernel Kii{t, t') given by (fT8|) . 
The Fourier transform of this kernel can be easily evaluated 

/oo -r 
dt e'^' K^it, t') = - cos{kt')e{a^ - k^) , (73) 
-oo 2 

with the step function 9{x) equal to 1 for x > and otherwise. Making use of this identity we 
perform Fourier transformation of both sides of (1351) and find that the function x(^) introduced 
in (17T1) vanishes for 6^ > B^ with B = aTi/2. This property should not be surprising since, in the 
Bethe Ansatz approach to the scaling function [UlTl], the function x(6') describes the distribution 
of hole^ which condense on the interval [— -B, B]. For —B < 9 < B the function x{9) is given by 

Ix {0) = l + m -^fdt (7h(t) - e*/^ Jo(2,t)) , (74) 



^We note that the absence of 0{p) term in the expansion of the scahng function, £2(3) = 0, is related to 
vanishing of 0{p^) term in the energy of nonrelativistic Fermi gas. 

^This function differs from the one introduced in [IT by normalization, x{(^) = "^jPhik) with 6 = fc7r/2. 
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where 9 = A;7r/2 and the notation was introduced for 



Let us rescale the integration variable, t t/2g, and apply (l39i) to eliminate ■y(t) = 7+(t) +7_(t) 
in favor of T±{t) 

(76) 

We already encountered similar integral in Sect. 3.2. The important difference is that the func- 
tions T±{t) are now defined for j ^ 0. This does not affect however the general scheme that 
we followed in Sect. 3.2 and calculation goes along the same lines as before (see Appendix B for 
details). It leads to the following relation in the scaling limit g —>■ oo and j/m = fixed (with 
9 = k'n/2) 

1 77 , ^ /•°° , cos(A;t)e*/2 1 / ,/o , , ? , A 
^^^^ = "2 + 8""'°'^^ + ^ X Ae-^ ^(^e*/^7h(t)-| Jo(2^t)J . (77) 

It is important to emphasize that this relation was derived in Appendix B under assumption that 

B <7ig. (78) 

Taking into account the identity 9 = kn /2, we find that this condition can be reformulated as a 
requirement that the interval [—a, a] on the real /c— axis in which the holes condense should be 
located inside the gap [—2g, 2g] corresponding to the magnon density at strong coupling described 
by the BES equation. 

Substitution of ([77]) into ([74]) yields 

16 f°° 1 + e* 2? f°° sinhft/2) 

X{9) =mcosh9 + — dt cos{kt)- ^7h(i) + — / dt cos{kt) M^gt) . (79) 

TT Jq 1 -|- e vr Jq cosh t 

The last term on the right-hand side of this relation is subleading in the scaling limit (see 
Eqs. ( IB. 151) and ( IB. 161) ) and can be neglected. Then, we replace 7h(^) by its Fourier integral 

2 f°° 1 + e* 

X{9) = m cosh 9 + — d9' x{9') / dt cos{29t/7i) cos {29't/7i) , (80) 

J^B Jo 1 + e 

integrate over t with a help of ( l68l) to find that x{^) = x(~^) satisfies the same integral equation 
( l66l) as the rapidity distribution in the 0(6) model! 

Finally, it follows from (l65l) and (!72|) that the density of particles in the 0(6) model is related 
to the scaling parameter j as 

P=7^ f d9x{9)=^im = L (81) 
271 J_B 2 

where in the last relation we applied fl24l) . 
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4.3 Energy of the ground state 

It remains to show that the scahng function e{g,j) is related to the energy of the ground state 
of the 0(6) model fl65|) through relation flTOj) . 

As follows from (l37j) . the scaling function admits the following representation 



di jf\2gt)-j'^>{2gt) 
t sinh(t/2) 



(0)/ 



(7h(t)-7h(0))+jei(^) 



^2) 



where we separated terms linear in j into the function ei{g) given by ( H9i) . As before, we eliminate 
'j±\t) in favor of r|^'*(t) with a help of (l39l) and replace 7h(i) by the Fourier integral (1721) 



e{9,j) = - [ dex{e)E{e)+je^{g), 

TT J-B 

where the notation was introduced for (with 6 = kir /2) 

.(«)^ff(i-co,H/.«,)[£^(r-,.,.rf,, 

-|^(^^°'<*)-^'w 



^3) 



M) 



We observe similarity of this integral with the one entering (176|) . Going through the same steps 
as before, we get at large g, 

TT) 

E{9) = — (cosh 9-1), (85) 

with m being the mass gap fl59l) . Plugging this expression into fl83l) and making use of fl56l) . we 
evaluate the scaling function as 



e(g^j) = — dex{0){coshe-l)+ji-l + m) = — d9xi0) coshe -J 

TT J-B TT J_B 



^6) 



in a perfect agreement with (ITOll and (l65l) . 

Thus, we demonstrated that, in the limit g ^ oo with j/m = fixed, the scaling function e{g,j) 
is related to the energy density in the ground state of the two-dimensional 0(6) sigma-model 
(l7Up and, therefore, it can be found from the exact solution of this model constructed in |2B]. We 
should keep in mind however that this result was based on the assumption (1781) which need to 
be checked. Let us first consider the region j <^ m in which case a = j /m + . . . and, therefore, 
B = an/ 2 automatically satisfies flTHl) . The relation flTHl) becomes nontrivial for j ^ m. In this 
region, the exact solution of the 0(6) model leads to 



eo(6) 



TT 



8 ln(j'/m) 



O 



In In (j'/m) 
In^ (j/m) 



^7) 



This expression is valid for j ^ m independently on whether j ^ g or j <^ g whereas the string 
theory consideration [2] suggests that the relation between e{g,j) and eo(6) should only hold for 
j <^ g. Remarkably enough, the relation ( 178|) leads to the condition for j on the gauge theory 
side compatible with j g. Indeed, for j ^ m we expect from the exact solution of the 0(6) 
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model [26l [271 EH] that the 5— parameter scales as i? ~ ln(j'/m) — clnln(j'/m) (with constant c 
that need to be calculated) leading to 



] r^mg 



where in the second relation we used the expression for the mass gap (l59ll . We conclude that the 
condition (1781) implies that j is exponentially small in g. We expect that the constant c should be 
equal to c = 3/4 leading to the condition j/g -C 1, in agreement with string theory prediction. 



5 Conclusions 

In this paper, we have studied anomalous dimensions of high-twist Wilson operators in a par- 
ticular limit [1] when their Lorentz spin grows exponentially with the twist. The anomalous 
dimensions have a nontrivial scaling behavior in this limit and conjectured integrability of the 
dilatation operator in A/" = 4 SYM theory leads to the FRS equation [H] for the corresponding 
scaling function e{g, j). It is believed that the function e{g, j) defined as a solution to this equation 
should interpolate between the weak coupling result based on explicit perturbative calculation 
in A/" = 4 SYM and the strong coupling result predicted within the gauge/string duality. 

We solved the FRS equation at strong coupling in the scaling limit, g ^ oo and j/m = fixed, 
and found that, quite remarkably, the anomalous dimensions of high-twist operators in four- 
dimensional conformal invariant A/" = 4 SYM theory are described by two-dimensional nonlinear 
0(6) model which is asymptotically free at short distances and develops a mass gap in the infrared. 
This result is in a perfect agreement with the proposal by Alday and Maldacena [2] who argued, 
based on the dual string consideration, that the scaling function e{g,j) should coincide with the 
energy density of the 0(6) model embedded into the AdSs x S^ model. 

In the Bethe Ansatz approach, the scaling function is determined by two sets of rapidities 
describing magnons and holes [1]. In the scaling limit, the rapidities condense on the real axis 
and their distribution densities satisfy the system of coupled integral equations [H]. We found 
that solutions to these equations at strong coupling depend on a 'hidden' nonperturbative scale 
m which determines the mass spectrum of hole excitations. The same scale has been previously 
identified in the strong coupling expansion of the cusp anomalous dimension [25]. We demon- 
strated by explicit calculation that, firstly, the scale m coincides with the exact mass gap of the 
two-dimensional 0(6) sigma- model and, secondly, the dynamics of holes at strong coupling is 
described by the thermodynamical Bethe Ansatz equations for the same model. 

At weak coupling, e{g,j) is a bi-analytical function of g and j [H]. At strong coupling, due to 
appearance of the mass gap, e{g,j) has a nontrivial 'phase diagram' in the {g, j)— plane. Namely, 
for m <^ g, the function e{g,j) has different behavior for j ^ m ^ g and m -C j -C [2]. 
In both regimes, e{g,j) coincides with the energy density in the ground state of 0(6) model: 
For j m g, the function e{g,j) admits an expansion in powers of j/m which reflects 
nontrivial properties of the 0(6) model in the infrared. For m <^ j g, in accordance with 
the asymptotic freedom in this model, e{g,j) is given by the perturbative series in 1/g with 
coefficients proportional to and enhanced by powers of \n{j/g). 

We would like to stress that the FRS equation predicts the scaling function e{g, j) for arbitrary 
g and j. The 0(6) sigma-model only describes the leading asymptotic behavior of e{g,j) in the 
limit g —>■ oo and j g. It would be interesting to investigate whether the relation with the 
0(6) model also exists at the level of sub leading corrections. Based on our experience with the 
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cusp anomalous dimension [22], |42l [35l [37l [29l [38l [39j, we expect that the transition from strong 
to weak couphng regime in e{g,j) should occur for g ^ 1. In this region, the hierarchy between 
the scales m and g disappears and the relation to the 0(6) sigma-model is lost. The crossover 
between strong and weak couplings deserves additional study. Finally, the scaling function €{g,j) 
has a nontrivial behavior [Ij at strong coupling in the region j > g (with N,L ^ oo), which 
is left beyond the scope of the present paper. As was already mentioned in the Introduction, 
the anomalous dimension in this region does not grow logarithmically with anymore but has 
instead a BMN like behavior ([3]). It would be interesting to reproduce this behavior from the 
FRS equation. 
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Appendix A: Derivation of the scaling function 



This appendix contains a derivation of the relation ( |371) . According to (|T5l) and ( 13T1) . the scaling 
function e{g,j) is given by 



e{g,j) = f{g,j) - f{g,o) = 16/ [71(^7,7) - 71(^7,0)] 



(A.l) 



where 71(5', j) enters into the Bessel series expansion of function '~f^{t;j), Eq. (l28l) . The function 
j-{t;j) satisfies the integral equations flM|) and satisfies Wronskian-like relations. 

To derive these relations we choose some reference j', multiply both sides of the two relations 
in fl34l) by the coefficients {2n — l)j2n-i{g, j') and {2n)'j2n{g, j'), respectively, and sum over n > 1. 
Then, we convert the sums into the functions 'j±(t;j') using the definition ( l28i) and subtract the 
second relation from the first one to obtain (after rescaling t t/2g) 



°° dt 
T 



^.{t-j)^^{t-j')-^+{t-j)^+{t-j') ^^{t-j)^+{t-j')+^+{t-j)^^{t-j') 

1 _ Q-t/2g + Qt/2g _l 

7i(^,j') + [{2n-i)h2n^i{gd)i2n-i{gd') - {'^n)h2n{gJh2n{g,j')]- 

n>l 



(A.2) 



The expression on the left-hand side is invariant under exchange j ^ j' and the same should be 
true on the right-hand side. Then, we replace hn by their definition fl5Ul) and get 



. .^ . .,^ 2 dt 7+(2^t;/)-7_(2(?t;/) 



J 



e'"^^{t-j)-'-U2gt) 



(A.3) 

where we indicated explicitly the dependence of the 7— functions on the scaling parameters. We 
note that for j' = the relation flA.Sp can be used to evaluate the scaling function flA.ip . We 
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take into account that 7h(i; j' = 0) = hn{g,j' = 0) = and '^±{t\j' = 0) = 7±'*(t) to obtain 



e{g,j) = ?>2g 



t e*-l 



e*/^h(t;j)-|)+i(l-^o(2^?t)) 



J 



(A.4) 



Here the integral involving the Bessel function can be further simplified by making use of the 
identity 1 - Jo(-2) = 2 Yl,n>i -hniz) leading to 



n>l 



dt^f{2gt)-^'^:^\2gt) 



(0), 



J2n{2gt) 



(A.5) 



-^93 



j^f\2gt){l-M2gt)) 



-^93 



dt (0) 



7r(2^t) 



where in the first relation we applied the integral equations for j = and in the second 
relation used Bessel series representation (128!) together with the orthogonality condition (132|) . 
Combining together the relations ( ]A.4I) and ( 1A.5I) we arrive at ( !37l) . 



Appendix B: Strong coupling expansion 



In this Appendix we present some details of the strong coupling expansion of the scaling function. 
In our analysis we encountered integrals like Eqs. (1531) and (17^ which involve the product of ratios 
of hyperbolic functions and linear combinations of T±{t). To begin with, we replace the former 
by their Fourier integrals with a help of identities (with > 0) 



cos(/ct) 
cos(A;t) 



cos\i{t/Ag) 
cosh(t/2^) 
sin\i{t/Ag) 
cosh(t/2^) 



\/2g I du cos('ut) 

J — oo 

\f2g / du sin(Mt) 



cosh(5f7r(M + A;)) 
cosh(25f7r(u + A;)) 

sinh(5f7r(u + A;)) 
cosh(2(77r(u + A;)) 



(B.l) 



valid for arbitrary real k and t. To verify these relations, one closes the integration contour into 
the lower or upper half-plane depending on the sign of t and picks up residues coming from the 
denominator. The calculation of 63 also involves the integrals 



t 



cosh(t/45() 
cosh(t/2^) 
^sinh(t/4^) 
'cosh(t/2^) 



-\f2q [ dusm(ut)—( 
7-00 du\ 



d f cosh.{giTu) 



V2g / du cos{ut) 



cosh.{2gTTu 
d / sinh(5f7ru) 
du \cosh(2(yf7rM) 



(B.2) 



In what follows we only calculate in detail ei since the analysis of es is essentially the same. 
Applying (]B.1|) for A; = we find from (153!) 



e^ig) = -2V2g' 



+ 



°° cosh (giTu) 

du — — 

cosh {2gTTu) 

00 



°° dt 

T 



cos (Mt)) [rL°Ht) + rf (t) 



(B.3) 



du 



sinh (giTu) 
cosh {2g7iu) 



dt 

— sm [ut) 



rW(t)-rV'^(t) 



(0), 
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Let us split the u— integrals into sum of two terms corresponding to < 1 and > 1. In 
the first one, we replace (1 — cos{ut))/t = dv sin(t>t) and sm{ut)/t = Jq dv cos(ft) and, then, 
evaluate the t— integral with a help of (15^ . The resulting expression for ei{g) is given by (l55ll . To 
find the large g expansion of the first term on the right-hand side of fl55|) we just do an opposite, 
that is rewrite the integral over — l<M<lasa difference of two integrals over —oo < u < oo 
and > 1 



sinh (gnu) 

duu — r 

I cosh [2g7TU) 



inu) 



°° , sinh 

du u — ^ '- 

cosh [2g'RU^ 

1 



^1 + Txg 
4^V!~ {-ngf 



-0(e 



duu 



-3g7r\ 



sinh (giTu) 
cosh {2g7Tu) 



(B.4) 



In the remaining two terms in fl55l) . we replace hyperbolic functions by their leading large g 
asymptotics, integrate over u by parts, take into account fl5^ to evaluate the boundary term and 
arrive at fl56|) with m given by 



m 



27/2 g J g-Trg 

TV \ ng 



1 

2 

X 



du e 



dt 



(B.5) 



cos(Mt) (rL°^(t) - rf (t)) + sin(Mt) {T^^\t) + rf (t) 



Integration over u leads to ( 157|) . Then, we change the integration variable in (157|) as t = (^vrr and 
take into account (!58|) to get 

27/2 



-Am , 



with Am given by a complicated twofold integral 

2^/2 /-"^ c/r e-^"/4r(3/4)r(l + ir/4) 



Am = ge~ 



■Kg 



71^ 



T + i 



r(3/4 + ir/4) 



-1 



u 



1 



1/4 



(B.6) 



(B.7) 



To find its large g asymptotics, we use Mellin-Barnes representation for e*'^^'^^ "^"^ and integrate 
over u to obtain 

-5+ioo 



Am 



-TTg 



TX 



— 5— ioo 



27rz^ r(j + 2) - 



^-zr)^'e-^'^/4r(l + ir/4) 
rT^ r(3/4 + ir/4) 



+ c.c. 



(B.8) 

At large g we deform the integration contour to the left and pick up the contribution of poles at 
negative j. The first two poles closest to the origin are located at j = —3/4 and j = —1 with 
the second one coming from integration at small r. Going through calculation of residues we get 



Am 



4 



{2.gr'l'''-^-{2.grH'l-^0[g-'^l') 



(B.9) 



Substitution of this relation into flB.6l) yields the desired expression for the mass gap fl59l) . 

The evaluation of the integral 1(6') defined in fl76|) goes as follows. We replace the factors in 
front of r±(t) functions on the right-hand side of ( !76l) by their Fourier integrals (see Eq. ( IB. II) ). 
split the M— integral into sum of two terms corresponding to -u^ < 1 (A+B) and u^ > \ (C) and, 
then, simplify the first term using ( H2l) . In this way, we find 



/(^)=JaW+/bW+/cW, 



(B.IO) 
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where 



9 

V2 



du 



cosh (qitu 



cosh {2g7iu + 29) ' 



(B.ii: 



/b(^) = 2V2g 



du 

1 ^0 



dt 



cosh(^7rM + ^) sinh(^7rM + ^) . , 
cos [2gut) H — — - sm [2gut) 



cosh {2gTTu + 26) 



cosh (25(7rM + 26) 



HO) 



2\f2 



du 



dt 



cosh {gnu + 6) 
cosh {2gnu + 26) 
sinh {gnu + 6*) 



cos(Mt) (^r_(t) - r+(t) 
sin(ut) fr_(t) + r+(t) 



cosh {2gnu + 26) 

Here Ia{0) and /b(^) originate from the first and the second terms on the right-hand side of f H2|) . 
respectively. 

To find asymptotic behavior of Ia{6) at large g we split the m— integral into the difference 
of two integrals over — oo < m < oo and u^ > 1. The former integral can be easily done with 
a help of identity (1B.1|) . Calculating the second integral as well as Ic{6)^ we replace the ratios 
of hyperbolic functions by their leading asymptotic behavior at large g assuming that {gnu ± 6) 
does not change sign inside the integration region. This assumption is justified provided that 
—ng < 6 < ng, or equivalently —2g < k < 2g. Since the integral (IB.lOp is defined on the support 
interval —B < 6 < B of the hole density x{0), the above mentioned condition is automatically 
satisfied for B < ng, or equivalently a < 2g. Assuming that this condition is fulfilled, we find 
after some algebra 

lA{e) + lc{0) = + ^m(j) cosh^ , (B.12) 

Z o 

where the parameter m{j) is given by 

= ^e--^-^ / due-^3u (B_13) 



n 



V2 



X 



cos(Mt) (r_(t) - r+(t) ) + sin(Mt) (r_(t) + r+(t) 



Here the dependence on j resides inside the functions T±{t) only. In the previous section, we 
encountered the same expression for j = and identified m(0) = m as the mass gap of the 0(6) 
sigma-model, Eq. (l59l) . We notice that the integral on the right-hand side of ( IB. 1311 scales as 
e~'^^ at large g. This allows us to neglect in the t— integral corrections proportional to j in the 
scaling limit g oo and j/m = fixed, that is to substitute the functions T±{t) by their values 
at j = 0, leading to = m(0) + . . .. From the point of view of the integral equation fH2l) this 
amounts to neglecting the second term on the right-hand side of P2|) as being proportional to j. 

To evaluate Ib{6), we rewrite the u— integral as the difference of two integrals over —oo < u < 
CO and > 1. The former integral can be easily evaluated with a help of identity fIB.ip while the 
latter one produces the contribution subleading in the scaling limit provided that —ng < 6 < ng 
and, therefore, can be safely neglected leading to 

e*/' W./2 .^ J 

'0 



Ib{0) 



dt 



cos{26t/n) 



e* -|- e" 



e*-l 



Jo{2gt) 



(B.14) 
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Combining together the relations ( IB.IOI) . (1B.12I) and (IB. 141) we arrive at (|77ll . 
Calculating (!79i) we encountered the following integral 

2j r , .o., / ^sinh(^/2) , y2j p , sinh(7r^M + ^) p , . , , 

/ dtcos{29tn) U^j^(2gt) = ^ du ^ , / rft Jot smut 

TT Jo cosht TT J_oo cosh(27r(7u + 26') Jo 

^ r ..^ sinh(7r^^ + 6^) - 1) ^ 

TT cosh(27r(7u + 26) v^^J^ ' ^ ' ^ 

where we rescaled the integration variable t —>■ t/2g, applied the second relation in fIB.ip and 
used the properties of Bessel functions. Expanding hyperbolic functions at large g we evaluate 
the integral as 

-^cosh^ / =-j(7-^/^e~-gcoshg + ... (B.16) 



vr Ji - 1 

In the limit g —>■ oo and j/m = fixed, it scales as and, therefore, provides a subleading 
contribution to the right-hand side of ( |79l) . 
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